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Abstract. Measurements of the pp — ► pK + A reaction at T p = 2.28 GeV have been carried out at COSY- 
TOF. In addition to the Ap FSI and TV* resonance excitation effects a pronounced narrow structure is 
observed in the Dalitz plot and in its projection on the p/l-invariant mass. The strongly asymmetric 
structure appears at the pp —¥ NK + £ threshold and is interpreted as UN cusp effect. The observed 
width of about 20 MeV/c 2 is substantially broader than anticipated from previous measurements as well 
as theoretical predictions. Angular distributions of this cusp structure are shown to be dissimilar to those 
in the residual pK + A channel, but similar to those observed in the pK + E° channel. 

PACS. 13.75.Cs - 13.75.Ev - 14.20.Jn - 14.20.Pt - 25.10.+S - 25.40.Ep 
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, 1 Introduction duces a non-analyticity in the spectra, a so-called cusp, 
£SJ ■ the strength of which is a measure of the corresponding 
t— I | The hyperon production in nucleon-nucleon collisions has transition matrix element (here EN — > AN) - e.g. see 
^ attracted interest in recent years for a number of reasons. Ref. [13] for the method to extract the s-wave 7T7t scat- 
First, it offers a valuable tool for the determination of tering lengths from K — > Sir. Signs of such a cusp in the 
the hyperon-nucleon final state interaction (FSI). Most di- pp — > pK + A reaction have been first observed in inclusive 
5h ■ rectly this effect is seen in the hyperon-nucleon invariant measurements with single-arm magnetic spectrometers at 
mass spectrum. For a recent determination of the ylp-FSI Saclay [T^] and COSY-HIRES [14J. A theoretical calcula- 
by the COSY-HIRES collaboration using a high-resolution tion of this cusp effect in proton-proton collisions has been 
magnetic spectrometer see Ref. [T|. Second, it gives access presented by Laget |15j . 

to the rare decay branches of N* resonances produced in Originally a peak structure at the position of the EN 

the NN collisions. For recent work on that see Refs. [3 threshold has been discovered in K~ absorption in deu- 

IH5IIS]. Third, it offers the chance to search for more teriumpES]. Subsequent bubble-chamber measurements [TH 

exotic objects like pentaquarks 0[EJ[H] or dibaryons [B [Jg] observed a pronounced peak at the E threshold with a 

HTJ. Among the latter an inevitable, though not very ex- width of about lfJ Me V/ c 2 and below. For a review and a 

otic candidate would be the EN system produced near discussion f possible dibaryon aspects see, e.g. Refs. [13 

threshold followed by a EN — > AN transition. Whenever [^rjir^nr^l 

an inelastic channel opens in a production reaction, it pro- ^ 
~ In exclusive and kmematically complete COSY-TOF 

a present address: Forschungszentrum Jiilich, D-52428 measurements indications of this cusp have been observed 

Jiilich, Germany at several energies [UEE ■ Since these measurements lack 

Correspondence to: H. Clement the necessary statistics for a reliable investigation of the 

email: clement@pit.physik.uni-tuebingen.de cusp effect, we use here the measurements at T p = 2.28 GeV 



2 



S. Abd El-Samad et al.: On the EN cusp in the pp — > pK + A reaction 



(p = 3.081 GeV/c) for a detailed investigation of this mat- 
ter. The primary purpose of this run, which comprises 
an order of magnitude higher statistics than the previ- 
ous TOF-measurements, was originally the pentaquark 
search [TJ. Angular distributions obtained from this run 
have been published already for both pK + A and pK + E° 
channels [2 J . The results presented in this paper are based 
on the thesis work of Refs. [231124] . where also details of 
experiment and analysis are found. 



2 Experiment 
2.1 Detector setup 

Since the experimental setup was discussed in detail al- 
ready in Refs. 0H1E], we give here only a short account. 
The measurements were carried out at the Jiilich Cooler 
Synchrotron COSY using the time-of-flight spectrometer 
TOF located at one of its external beam lines. The TOF 
spectrometer is a modular detector setup, which can be 
adapted to the specific requirements of an experiment. 
Here it was used in its standard version for hyperon pro- 
duction, see e.g. Figs. 1 and 2 of Ref. [I]. At the entrance 
to the detector system the beam - focused to a diame- 
ter smaller than 2 mm - hits the thin- walled LH2 target, 
which has a length of 4 mm, a diameter of 6 mm and 0.9 
pm thick hostaphan foils as entrance and exit windows. 
At a distance of 22 mm downstream of the target the 
two layers of the start detector (each consisting of 1 mm 
thick scintillators cut into 12 wedge-shaped sectors) were 
placed. A silicon microstrip detector as well as two fiber 
hodoscopes were installed at distances 30, 100 and 200 mm 
from the target. These three tracking detectors provide 
the position information of the traversing charged parti- 
cles, whereas the start detector supplies the start time 
signals for the time-of-flight (TOF) measurements. After 
a flight path of about 3 m through the evacuated vessel 
the charged particles are detected in the highly segmented 
stop detector system consisting of the triple-layered quirl 
and ring detectors as well as a single-layered 96-fold seg- 
mented barrel detector. 



2.2 Particle identification and event reconstruction 

In the experiment the trigger suitable for the selection of 
hyperon production events required two hits in the start 
detector and four hits in the stop detectors. This multiplic- 
ity jump from two to four specifically selects the produc- 
tion of neutral hyperons, which decay into charged prod- 
ucts like the A decay process A — > pir~ , which happens 
with a branching fraction of 64%. Tracks of charged par- 
ticles are reconstructed from straight-line fits to the hit 
detector elements in start, fiber and stop detectors. They 
are accepted as good tracks of primary particles (p, K + ), 
if they originate in the target. The secondary particles re- 
sulting from hyperon decay (p, it~ ) form a V^-shaped track 
pair originating from a secondary vertex downstream of 



the microstrip detector. Primary vertices, which are lo- 
cated within the target volume, were reconstructed with 
an accuracy of o~ XjV = 0.25 mm and a z = 0.7 mm. The 
secondary vertex from the A decaying downstream the 
microstrip detector and upstream the first hodoscope was 
reconstructed with accuracy of a x ,y = 1.5 mm and o~ z = 
4.5 mm. 

In addition to the tracking information we use the TOF 
information of all four ejectiles of an event to determine 
their four-momentum vectors. For the particle identifica- 
tion a kinematic fit is applied, where all permutations of 
particle assignments are considered and the one with the 
best x 2 is selected as the correct one. Since the kinematics 
of the light (7r, K) and heavy (p) emitted particles is quite 
different, finding of the correct particle assignments actu- 
ally is rather clear-cut. According to detailed Monte Carlo 
(MC) simulations the misidentification rate was in the or- 
der of a few percent only - in agreement with previous 
TOF results 00]. 

Fig. 1 depicts the two-dimensional plot of the invariant 
mass M p7T - of the A decay particles versus the primary 
particle pK + -missing mass spectrum for MC simulation 
(top) and data (bottom). As can be seen the resolution in 
invariant and missing mass are comparable, and the MC 
simulation reproduces very well the experimental situa- 
tion. The dashed circles in Fig. 1 indicate the region of 
events accepted for the subsequent analysis steps. 

The kinematic fit was fivefold overconstrained, where 
the fifth overconstraint originates from the condition that 
the invariant mass of the four-momenta of the decay par- 
ticles has to be identical to the A mass. A total of 30,000 
events passed the criterion of prob(x 2 ) > 10% in the kine- 
matic fit for being accepted as a proper event [23]. 

The primary particle pK + missing mass spectrum of 
the events finally selected by the % 2 criterion is shown in 
Fig. 2 before the kinematic fit. The A peak appears with 
essentially no background and is in very good agreement 
with the Monte Carlo (MC) simulations of the detector 
performance. The slight mass shift of the peak between 
data and MC is due to imperfections in the energy cali- 
bration. 

The MC simulations also have been used for efficiency 
and acceptance corrections of the data. We used for these 
simulations as input a model-description for the reaction 
of interest, which provides a good description of the exper- 
imental differential distributions - see discussion in section 
4. We note, however, that since the TOF detector covers 
nearly the full phase space of the reaction - see the Dalitz 
plots in Fig. 3 - the differences in the corrections between 
model-based and pure phase MC simulations are only mi- 
nor. 

By use of the kinematic fit the mass resolution in 
the spectra of the invariant masses Mr a, M p k and M p a 
(Fig. 4) improves from about 30 MeV/c 2 to 6 MeV/c 2 
(FWHM). 

The absolute cross section is obtained by relative nor- 
malization to literature data [25 , 26J for the elastic pp 
channel, which was measured in parallel during the ex- 
periment. 
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Fig. 1. Plot of the invariant mass M p7V - := Mdecay 
decay particles versus the primary particle pA" + -missing mass 
MMpK for MC simulation (top) and data (bottom). The 
dashed circles indicate the range accepted for the subsequent 
kinematic fit. The scale of the z-axis is in arbitray units. 
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Fig. 2. p_ft' + -missing mass spectrum of events, which were 
finally selected as good candidates for the pK + A channel by 
the x 2 criterion. The shaded area gives the MC simulation, 
the vertical middle line the position of the A mass and the 
two outer vertical lines the range accepted for the subsequent 
kinematic fit. 



Single differential cross sections are shown in Figs. 4 - 
6. Fig. 3 exhibits the acceptance and efficiency corrected 
Dalitz plot of the three-body exit channel pK + A. It is by 
no means homogeneous, i.e. phase space- like. The inten- 
sity in the Dalitz plot peaks at the left side of the short di- 
agonal corresponding to high pK + -invariant masses M pK + . 
In addition we see two vertical narrow structures corre- 
sponding to the pyl-invariant masses M p a at threshold and 
at M pA w 4.54 GeV 2 /c 4 (see vertical arrows in Fig. 3), i.e. 

n p + ms- The first structure may 



M pA ~ 2.13 GeV/c 2 « r, 
be related to the Ap FSI, whereas the latter one is in the 
region of the Up production threshold. 

In Fig. 4 we show the spectra of the three invariant 
mass systems Mka, MpK and M p a The data of all three 
strongly deviate from phase space, which is indicated by 
the shaded area in the plots. Among them the M p a spec- 
trum appears particularly interesting, since it exhibits two 
narrow structures. They may be connected with FSI and 
£"N cusp as mentioned above and will be discussed in the 
following. 



3 Results 



4 Discussion 



For the total cross section of the pp — > pK + A reaction at 
2.28 GeV we obtain (21.2 ±0.2 ±2.0) ^b, where the errors 
are statistical and systematic respectively. Within uncer- 
tainties this value agrees with that obtained in Ref. [2] • As 
in the previous TOF results [2 ,4,5,6J for the pp — > pK + A 
reaction the systematic uncertainty of about 10% is the 
by far dominant uncertainty. It originates from the uncer- 
tainties in the luminosity determination, reconstruction 
efficiency and acceptance correction. 



As shown in previous works PJIUEJE], the gross features 
of the pp — > pK + A reaction may be well described by N* 
resonance production with subsequent N* — > K + A de- 
cay and by inclusion of the Ap FSI. The A^* resonances, 
which play a role here, are ^(1650)1/2", A^(1710)l/2+ 
and AT(1720)3/2+. For the description of the data we as- 
sume excitation of these N* resonances via t-channel me- 
son exchange along the prescription given in Refs. [4"ll2"3"]. 
The £ cusp is treated in a simplified manner just as a nar- 
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Fig. 3. Dalitz plot of M^ +A versus M pA for the pp — s- pK + A 
reaction at T p = 2.28 GeV, for MC simulation (top) and data 
(bottom). The two vertical narrow structures (see arrows) at 
the pA threshold and at M pA = 4.54 GeV 2 /c 4 are due to the 
pA FSI and the EN cusp, respectively. The scale of the z-axis 
is in arbitray units. 



row Ap resonance and the Ap FSI is taken into account in 
the factorization approximation of Ref. pQ. 

For the model fit to the data we allow mass and width 
of the N* resonances to vary within the boundaries given 
in PDG [57|- I n addition relative phase and strength pa- 
rameter for each of the resonances have been fitted. For a 
quantitative description of the data in the FSI region we 
also need a readjustment of the FSI parameters resulting 




M KA [GeV/c 2 ] 




l\^ K [GeV/c*] 




IV^JGeV/c] 



Fig. 4. Differential distributions of the invariant-mass sys- 
tems M K + A (top), M pK + (middle) and M pA (bottom). The 
shaded areas indicate phase-space distributions, the dotted 
lines the contributions of iV* resonances, the dashed line the 
E cusp effect and the solid line the full MC simulation. 
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Table 1. Results from the fit to the data for N* resonances. 



resonance 


mjv 


r N * 




a 




(MeV) 


(MeV) 


(deg) 


/Jo 


iV(1650)l/2- 


1653 


168 


29 


3.5(4) 


AT(1710)l/2+ 


1712 


99 


16 


3.5(4) 


iV(1720)3/2+ 


1731 


383 





12.8(13) 



in a Ap = -2.2 fm and r\ v = 1.4 fm. These values are within 
la and 2a, respectively, of the HIRES result pQ. 

The result of a fit to the data is shown in Fig. 3, top, 
as Dalitz plot and in Fig. 4 by dotted (TV* resonances), 
dashed {£ cusp) and solid (full calculation) lines. The re- 
sulting values for mass, widths and relative phase of the 
N* resonances and S cusp are given in Table 1 in ad- 
dition to their total cross section contributions. The data 
are reasonably well described by this fit. According to this 
analysis the enhancement (relative to phase space) in the 
M p A spectrum at the pA threshold is predominantly due to 
the Ap FSI, whereas the enhancement at high pA masses 
arises from N* excitation, in particular from the excita- 
tion of 7V(1650)l/2~ - see Fig. 4. A much more sophis- 
ticated Dalitz plot analysis in the framework of Ref. [28J 
taking into account all COSY-TOF data on this reaction 
is in progress. Hence we will not discuss the contribution 
of N* resonances here further, but rather concentrate on 
the discussion of the cusp effect. The primary goal of the 
fit here is just to have a reasonable description of the data 
for the purpose of a reliable acceptance and efficiency cor- 
rection by MC simulations - as discussed above in section 
2 - and to also have some reliable estimate of the physi- 
cal background (due to N* resonances) underneath the £ 
cusp. Also, as we see from the fit in Fig. 4, bottom, the 
assumption of a Breit-Wigner distribution for the £ cusp 
is not a good description of this phenomenon. 

In fact, since the cusp is a threshold phenomenon we 
do not expect a symmetric Breit-Wigner distribution, but 
an asymmetric energy dependence in form of a Flatte dis- 
tribution \2~§\ 



da/dM Ap 
with 



r A P I \m% - m\ - im R (r Ap + r Sp )\* 



r Ap = gApQAp and r Sp = gs P qsp, 

where gi and qi are coupling constants and cm momenta, 
respectively, in the corresponding two-body subsystems. 
We have 

•J ( m |p-( m ^+ m p) 2 )(™i p -( m P- m ^) 2 ) 
IZp ~ 2^71 



and 



y j / ((m s +rn p ) 2 -rn 2 Sp )(rn 2 IJp -(m p -m I :) 2 ) 



above and below threshold, respectively. 

If 9en *C 9Ap then the Flatte distribution approaches 
a symmetric distribution. If gsN ^ 9Ap we have a very 
asymmetric distribution with a trailing slope at energies 
below the cusp and a rapid decline beyond the cusp. This 
situation is just opposite to what we observe in our data. 
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Fig. 5. Differential distribution of the invariant-mass system 
M p a- The shaded area represents the phase-space distribution, 
the solid curve shows the shape of the cusp as obtained from a 
coupled-channel treatment [30] . averaged over the experimen- 
tal resolution and fitted in height to the data at the low-energy 
side of the cusp. The cusp is assumed to sit upon a smooth 
background represented by the dashed line. Top: data of this 
work taken at 3.081 GeV/c. Middle: data taken with the new 
straw tracker at 2.95 GeV/c, normalized arbitrarily |24| . Bot- 
tom: comparison of both data sets in the cusp region. 
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Also, a large g% n means a strong EN FSI, which in 
turn causes a strong low-mass enhancement in the EN 
invariant-mass spectrum as well as a steep increase of the 
total cross section near threshold of the pp — > ENK reac- 
tions. Explicit measurements of these E production chan- 
nels [3J[3J provide no evidence for that. Since we have a 
sizeable Ap FSI, we hence expect g^N "C gAp, & more 
or less symmetric distribution around the cusp. This is 
borne out also in a recent theoretical treatment of A and 
E production in AW collisions in the framework of the 
coupled channel effective range method [30] ■ The resulting 
cusp distribution folded with the appropriate experimen- 
tal resolution is shown in Fig. 5 as solid curve. Similar to 
the prediction of Laget [IS] we have the situation that the 
low-energy side agrees very well with the data, whereas 
the data on the high-energy side fall off much less steep 
than predicted. 

In order to check whether the observed larger width 
is due to our experimental resolution, we compare in Fig. 
5 our M p A data from the high-statistics run with a more 
recent COSY-TOF high-resolution measurement at 2.95 
GeV/c [23], see Fig. 5, middle. This measurement was per- 
formed utilizing the new straw-tracker at COSY-TOF pro- 
viding an invariant mass resolution of 2.6 MeV/c 2 FWHM, 
i.e. more than twice better than in this work. Both data 
sets are compared in Fig. 5, bottom, depicting the cusp 
region in enlargement. Aside from a possible slight shift 
of about 1-2 MeV/c 2 , which is within the uncertainty of 
mass calibrations, both data sets coincide - in particular, 
if the different resolutions are taken into account. We see 
that both data sets exhibit the cusp in compatible shape 
and also size (relative to the background). 

Admittedly, there appears to be possibly a slight dif- 
ference. Whereas the high-statistics data exhibit a gentle 
decline at the high-energy side, the high-resolution data 
show the indication of a roughly 5 MeV broad bump upon 
the declining slope near 2.146 GeV/c 2 . Since this bump 
effect may be less than 3<t - depending on the assumption 
of background - and hence not statistically significant, we 
do not want to speculate about its nature at this point. 
However, we would like to mention that already kaonic 
deuterium data suggested a two-bump scenario. In partic- 
ular the bubble-chamber data with the highest statistics 
and a quoted mass resolution of cr~ 1.0-2.6 MeV/c 2 [T7] 
have been fitted by two symmetric Breit-Wigner functions 
with mi = 2128.7 ± 0.2 MeV/c 2 , A = 7.0 ± 0.6 MeV/c 2 
and m 2 = 2138.8 ± 0.2 MeV/c 2 , r 2 = 9.1 ± 2.4 MeV/c 2 
suggesting two nearby resonance states separated by 10 
MeV/c 2 . The position of this second bump (shoulder) is, 
however, not compatible with the position of the uninci- 
sive bump at 2.146 GeV/c 2 in Fig. 5, middle and bottom. 
From this we conclude that at present there is no statis- 
tically solid evidence for a second bump beyond the EN 
cusp, however, there is a solid evidence for a surplus of 
cross section right beyond the cusp position, which so far 
is not understood theoretically. 

The situation might change and the physical relevance 
of the second bump structure seen in Fig. 5, middle, might 
have to be re-discussed, if it should turn out that the cusp 



effect has indeed an energy dependence at the high-energy 
side as given by Refs. |1 51130] . Having a width of only a 
few MeV, such a bump appears to be smeared out in our 
measurement with the coarser energy resolution shown in 
Fig. 5, top, producing there just a shoulder. Also, the dif- 
ferences between the two-resonance scenario of Tan et al. 
[1 TJ and our results for cusp and second bump could be 
reconciled, if we allow to shift their resonance masses by 
3 MeV/c 2 , so that their first resonance coincides with the 
value for the EN threshold - as we find it for the position 
of the cusp in our data. That way the value for the second 
resonance mass in Ref. [T7] would move to 2.142 GeV/c 2 
and be no longer in serious disagreement with the posi- 
tion of the second bump in our high-resolution spectrum. 
For a new reinvestigation of this two-resonance scenario 
including the data of this work we refer to Ref. [33J . 

The different behavior at the high-mass end of the two 
data sets displayed in Fig. 5 is just due to the different 
beam energies: For the 2.95 GeV/c data the kinematic 
phase-space limit is 2.258 GeV/c 2 , for the 3.08 GeV/c data 
the corresponding value is 2.299 GeV/c 2 . In addition the 
influence of the broad N* resonances, which determine the 
continuum below the cusp, depends on the beam energy 
- as demonstrated in Ref. [3]. 

As expected the cusp is located at m p + m^o = 2131 
MeV/c 2 . This is consistent with the assumption that the 
observed structure corresponds to the production of E 
right at its threshold. In principle there should be two EN 
cusps, since m p +m s a = m n +mx:+ + 2 MeV/c 2 . However, 
due to the finite invariant-mass resolution of the data we 
are unable to separate these. Since the excitation of N* 
resonances gives a flat energy dependence in the M p a dis- 
tribution, see Fig. 4, we assume the physical background 
not belonging to the cusp scenario to be represented by 
the dashed curve drawn in Fig. 5. Associating that way 
the surplus of cross section above the dashed line with 
the cusp scenario, i.e. neglecting possible interference ef- 
fects, we obtain a cusp cross section, which corresponds 
to roughly 5% of the total pp — > pK + A cross section. 

The cusp scenario has a variety of consequences, which 
can be tested experimentally by the cusp angular distribu- 
tions. Such angular distributions are shown in Fig. 6. They 
have been obtained by subtraction of the background as 
shown in Fig. 5, however, now for each angular bin in- 
dividually. I.e., we split the double-differential cross sec- 
tion d 2 a / dMpAdcosO into eight equidistant angular bins 
of dcosO, fit then in each of these eight spectra a second- 
order polynomial to the data in the regions outside the 
cusp and integrate the strength in the cusp region above 
this polynomial line. This procedure of subtracting the 
background for each angular bin individually has been ap- 
plied successfully already in previous TOF work, see e.g. 
Refs. [H[3], for a detailed presentation of this method see 
Ref. [3T]. We note that a sideband background subtrac- 
tion [23J leads so similar results, however, it is in general 
not as reliable. 

Since the cusp effect is small compared to the back- 
ground originating from N* production, the cusp angular 
distributions have now substantial uncertainties, both sta- 
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Fig. 6. Distribution of the EN cusp over the K + angle in 
the center-of-mass system @2+ = 180° — 0^™ (top) and over 
the proton angle in the pA subsystem (Jackson frame) 0p A 
(bottom). The plotted error bars include both statistical and 
systematic uncertainties. The shaded areas denote phase-space 
distributions. The solid curve in the top figure represents a 
Legendre fit. 

tistical and in particular systematic. The systematic errors 
have been studied by varying the background description 
with the second-order polynomial. As a result we find that 
the systematic uncertainty is up to four times as big as 
the statistical uncertainty. The latter is readily estimated 
from the fact that the cusp is just 1/20 of the total cross 
section. With a total of 30000 events and eight angular 
bins we have on average only somewhat more than 100 
events left in the cusp of an angular bin spectrum, i.e. the 
statistical uncertainty is already in the 10% region. 

The condition of the E being produced right at thresh- 
old means that the E and N are in relative s wave with 
the consequence that the spin-parity of the EN system 
right at threshold must be J p = + or 1 + . Hence, the 
pA system resulting from the EN system can only be in 



relative s- or d-waves. The observed proton angular dis- 
tribution in the pA subsystem (Jackson frame) within the 
EN cusp is compatible with s— wave phase-space, see Fig. 
6, bottom. A dominant d-wave contribution would lead 
to a strongly anisotropic distribution. We note that the 
observed distribution is very different from the situation 
in the residual pK + A channel, where the corresponding 
Jackson frame distribution exhibits strong p-wave contri- 
butions - see Fig. 7 in Ref. [2J. 

The angular distribution of the K + , shown in Fig. 6, 
top, is compatible with dominantly s- waves relative to the 
EN system, which is observed in our case as pA system in 
the EN cusp. The p-wave component could be obtained 
from the Legendre fit to the data according to the ansatz 

da/dcosO^ 11 ~ a + a 2 (3 cos 2 6>™ 1 - l)/2, 

where the parameters ao = (0.50 ± 0.10) fib and a 2 = 
(0.28 ± 0.09) /ib. This result is different from the situa- 
tion in the residual pK + A channel, however, very similar 
to that observed in the pK + E° channel at the same in- 
cident energy, see Fig. 10 and Table 5 in Ref. [2J, and it 
is qualitatively similar to the K° cm angular distribution 
measured in pK°E + channel at 2.95 GeV/c |32| . 

As mentioned in the introduction the only other data 
on the EN cusp in pp induced K production originate 
from inclusive single- arm magnetic spectrometer measure- 
ments at Saclay [12J and COSY (HIRES collaboration 
|10|). In both cases essentially only a sharp increase in 
the cross section with a slight indication of a bump is seen 
in the K + missing mass spectrum at the Ep threshold. 
Since in these inclusive measurements E production can- 
not be separated from A production, the trailing slope of 
the EN cusp is not observed. If we compare these miss- 
ing mass spectra to the M pj \ spectrum in Fig. 5 - though 
the latter contains the integration over all K + angles - we 
see that at least qualitatively the shape of the single-arm 
spectra is very close to that of the M p a spectrum up to 
the maximum of the EN cusp. This is true also for the 
height of the cusp relative to the pK + A continuum left 
from the cusp. 

In Ref. [T2] the cusp has been extracted from data at 
T p = 2.3 GeV - which is close to our energy - by subtracting 
phase space distributions for A and E production - see 
Fig. 5 in Ref. [12J. As a result of this very crude treatment 
they obtain a bump at the cusp position, which indicates a 
width of roughly 10 MeV/c 2 - at variance with our findings 
of a much broader structure. 

Also, the width of 20 MeV/c 2 observed now in ex- 
clusive measurements is much larger than the value of 3 
MeV/c 2 assumed in Ref. [10J for the width of the EN 
cusp. Since the width of the EN cusp affects sensitively 
the value for the total cross section of the pp — > nK + E + 
cross section extracted in Ref. [TU] from the inclusive K + 
missing mass spectrum, the value of that cross section re- 
duces substantially when accounting for the much larger 
width of the cusp extreacted here. 

The calculations of Laget, which account for both A 
and E production as well as the EN cusp effect, are pub- 
lished only for specific K + scattering angles. This leaves 
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us merely with the possibility of a very qualitative com- 
parison. In Refs. [12,15] the calculations are shown for a 
K + scattering angle of 10° at T p = 2.3 GeV. Since the 
latter is close to the incident energy of this work, we may 
compare that calculation directly to the M p a spectrum in 
Fig. 5. From visual inspection we see that Laget's calcu- 
lation qualitatively gives the right order of magnitude for 
this effect, however, the calculated cusp effect appears to 
be substantially narrower than observed in our data — as 
already mentioned above. 

The cusp possibly may shed also new light onto the 
question about the ZW-FSI. From differential and total 
cross section measurements of the pp — > pK + E reac- 
tion it has been concluded that there is no sizeable E°p- 
FSI [2]. The presence of the pronounced cusp observed 
in pp — > pKA needs a different interpretation, since such 
a structure points to a very strong EN — > AN transi- 
tion. Therefore even at the threshold the UN interaction 
is strongly inelastic, which might well be the reason of 
diminishing any strong final state distortion in the exci- 
tation function for pp — > NKE. It is therefore particu- 
larly interesting to look into data on the pp — > pK°E + 
reaction, where the EN final state is purely isospin 3/2 
and therefore does not couple to the AN channel. In fact, 
differential and total cross section data for this channel 
have been taken recently at TOF [3J. They show no sign 
of any significant FSI effects. Following the argumenta- 
tion above this means that there is no sizeable FSI in the 
1=3/2 channel - different to the situation in the 1=1/2 
channel. A thorough theoretical explanation of these find- 
ings is highly desirable. For recent data on that channel 
see Ref. [3] and references therein. 



5 Summary 



This work presents the first exclusive and kinematically 
complete measurements of the EN cusp effect in the pp — > 
pK + A reaction and establishes this phenomenon in proton 
induced A production for the first time. The data exhibit 
a pronounced asymmetric shape of the cusp, which is gen- 
tly declining at its high-energy side. This is opposite to 
what is expected from a Flatte distribution. Due to the 
gentle fall-off at its high-energy side the cusp appears to 
be significantly broader than anticipated from theoretical 
predictions. Whether this is indicative of a narrow res- 
onance above the cusp energy - as speculated in earlier 
measurements of kaonic deuterium and as also possibly 
suggested by the high-resolution data presented here, can 
not be decided at the present stage. 

The measured angular distributions of the cusp point 
to s- waves between the kaon and the pE system, as well as 
dominantly s-waves between the subsequently emerging A 
and proton. Detailed theoretical calculations for this cusp 
effect would be highly welcome in view of these new higher 
statistics measurements. 
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